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Abstract 

The Muskat problem models the evolution of the interface between two different fluids 
in porous media. The Rayleigh- Taylor condition is natural to reach linear stability of the 
Muskat problem. We show that the Rayleigh- Taylor condition may hold initially but 
break down in finite time. As a consequence of the method used, we prove the existence 
of water waves turning. 

1 Introduction 

The Muskat problem [26] models the evolution of an interface between two fluids of different 
characteristics in porous media by means of Darcy's law: 

^u = -Vp-{0,gp), (1) 

K 

where {x,t) G x M"*", u = {ui{x,t),U2{x,t)) is the incompressible velocity (i.e. V • u = 0), 
p = p{x, t) is the pressure, /i(x, t) is the dynamic viscosity, k is the permeability of the 
isotropic medium, p = p{x,t) is the liquid density, and g is the acceleration due to gravity. 
More precisely, the interface separates the domains fi^ and ^l'^ defined by 



{p,p){xi,X2,t) 



and p} , p? , p^ , p^ are constants. This physical situation is also related to the evolution of two 
fluids of different characteristics in a Hele-Shaw cell |22j . due to the fact that the laws which 
model both phenomena are mathematically analogous |31| . 

This paper is concerned with the case p^ = p^ which provides weak solutions of the 
following transport equation 



pt + u- V p = Q, 
po = p(2;,0), x € 



(2) 
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where initially the scalar pQ is given by 

n -n(r r n^ - I ^'^^^ = ^^^^ > /o(^l)} 

Let the free boundary be parametrized by 

dn^{t) = {z{a,t) = (zi(a,t),Z2(a,t)) : a S M} 

where 

z{a, t) — {a, 0) 

is 27r-periodic in the space parameter a or, an open contour vanishing at infinity 

lim {z{a,t) - (a,0)) = 

a— >±oo 

with initial data z{a,0) = zo(a) = (a,/o(a)). From Darcy's law, we find that the vorticity 
is concentrated on the free boundary z{a,t), and is given by a Dirac distribution as follows: 

V"*" • u{x, t) = uj{a, t)6{x — z{a, t)), 

with uj{a,t) representing the vorticity strength i.e. • ti is a measure defined by 

<V'^-u,r]>= / uj{a,t)ri{z{a,t))da, 



with rj{x) a test function. 

Then z{a, t) evolves with an incompressible velocity field coming from the Biot-Savart 
law: 

u{x, t) = V-^A^^V-^ • u{x, t). 

As (x, t) approaches a point z{a, t) on the contour the velocity u agrees, modulo tangential 
terms, with the Birkhoff-Rott integral: 

BRiz,u)ia,t) = ^PV [ n4^4^4|S'^(^'*)^^- 
27r J \z{a,t) - z{(3,t)\'^ 

This yields an appropriate contour dynamics system: 

zt{a, t) = BR{z, Lo){a, t) + c{a, t)daz{a, t), (4) 

where the term c represents the change of parametrization and does not modify the geometric 
evolution of the curve |24j . 

The well-posedness is not guaranteed in general, in fact such a result turns out to be false 
for some initial data. Rayleigh [30] and Saffman- Taylor [31] gave a condition that must be 
satisfied for the linearized model in order to have a solution locally in time, namely that the 
normal component of the pressure gradient jump at the interface has to have a distinguished 
sign. This is known as the Rayleigh- Taylor condition: 

a{a,t) = -{Vp^{z{a,t),t) - Vp^(z{a,t),t)) ■ d^z{a,t) > 0, 
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where Vp' {z{a,t),t) denotes the hmit gradient of the pressure obtained approaching the 
boundary in the normal direction inside i}^{t). We call a{a,t) the Rayleigh- Taylor of the 
solution z{a, t). 

Understanding the problem as weak solutions of ([T]l2]) plus the incompressibility of the 
velocity, we find that the continuity of the pressure {p^{z{a,t),t) = p^{z{a,t),t)) follows as 
a mathematical consequence, making unnecessary to impose it as a physical assumption (for 
more details see [13] and [H]). For the surface tension case, there is a jump discontinuity of 
the pressure across the interface which is modeled to be equal to the local curvature times 
the surface tension coefficient: 

p^{z{a,t),t) — p^{z{a,t),t)) = TK{a,t). 

This is known as the Laplace- Young condition, which makes the initial value problem more 
regular. Then there are no instabilities [18] but fingering phenomena arise [29^ I19|. 

By means of Darcy's law, we can find the following formula for the difference of the 
gradients of the pressure in the normal direction and the strength of the vorticity: 

a{a,t) = {p^ — p^)daZi{a,t) 

uj{a,t) = -{p"^ - p^)daZ2{a,t). (5) 

Above g is taken equal to 1 for the sake of simplicity. 

Then, if we choose an appropriate term c in equation ^ (see section [2] below), the 
dynamics of the interface satisfies 

zt{a,t) = ^^PV [ - ^liM) _ d^z{P,t))dp. (6) 

ivr J \z[a,t) - z[P,t)\^ 

A wise choice of parametrization of the curve is to have daZi{a,t) = 1 (for more details 
see |13]). This yields the denser fluid below the less dense fluid if p^ > p^ and therefore the 
Rayleigh- Taylor condition holds as long as the interface is a graph. This fact has been used 
in |13j to show local existence in the stable case (p^ > p^), together with ill-posedness in the 
unstable situation (p^ < P^)- Local existence for the general case (/U^ 7^ /i^) is shown in |11] . 
which was also treated in \3A\ [1] . 

^From ([6|) it is easy to find the evolution equation for the graph: 

f,{a,t) = ^^^PV [ ^.^Ju^^l f(o ...2 {daf{a,t)-daf{f3,t))dl3, 

2vr J^{a- + [f[a,t) - f[/3,t)y (7) 

/(a,0) = /o(a). 

The above equation can be linearized around the flat solution to find the following nonlocal 
partial differential equation 

ft{x,t) = -^^^Af{x,t), 
/(x,0) = /o(x), xGM, 
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where the operator A is the square root of the Laplacian. This Hnearization shows the 
parabohc character of the system. 

Furthermore the stable system gives a maximum principle ||/||j;^oo (t) < ||/||j;^oo (0) 
decay rates are obtained for the periodic case: 

II/IIl-W < ll/ollL-e-^*, 
and also for the case on the real line (flat at infinity): 

" fo||L°° 



l + Ct 



There are several results on global existence for small initial data (small compared to 1 
in several norms more regular than Lipschitz [9l [35} [321 [T3l I19j ) taking advantage of the 
parabolic character of the equation for small initial data. In [8] it is shown in the stable case 
that global existence for solutions holds if the first derivative of the initial data is smaller 
than an explicitly computable constant greater than 1/5. Furthermore, if ||/o||loo < oo and 
IIc^x/oIIl"" < 1) then there exists a global-in-time solution that satisfies 

f{x,t) G C([0,r] X M) nL°°([0,r];H^^'°°(M)), 

for each T > 0. In particular / is Lipschitz continuous. 
Moreover, equation ^ yields an decay: 



,Xt) + — [ ds [da [ dx In (l + 

27r Jq Jm V 



7;Z / / / ^ I ^TT^;^ I I = II-/OIIL25 

which does not imply, for large initial data, a gain of derivatives in the system (see [8]). We 
will see below that the solutions to the Muskat problem with initial data in H'^ become real 
analytic immediately despite the weakness of the above decay formula. 
The main result we present here is: 

Theorem 1.1 There exists a nonempty open set of initial data in iJ^ with Ray leigh- Taylor 
strictly positive a > such that in finite time the Rayleigh- Taylor a{a,t) of the solution of 
([6]) is strictly negative for all a in a nonempty open interval. 

The geometry of this family of initial data is far from trivial: numerical simulations performed 
in |16j show that there exist initial data with large steepness for which a regularizing effect 
appears. In fact, as will be explained in Section [21 the first evidence of a change of sign in 
the Rayleigh- Taylor has been experimentally found in a model with two interfaces. 
We proceed as follows: 

First, in section 3, we assume initial conditions at time t = to that satisfy the Rayleigh- 
Taylor (a > 0) and the arc-chord condition, and for which the boundary z initially belongs to 
H^. Let Ci be the constant in the arc-chord condition, let C2 be an upper bound for the 
norm of the initial data and let C3 be a lower bound for a. Then there exists > to, with ti 
depending only on Ci, C2, C3, such that the Muskat problem has a solution for time t € [to, ii], 
satisfying also the arc-chord and Rayleigh- Taylor conditions. Moreover, for to < ^ ^ ^i, the 
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solution z{a,t) is real analytic in a strip S{t) = {a + : \C\ < c{t — to)}, where c depends 
only on Ci, 02,03. 

Our goal in section 4 is to show that the region of analyticity does not collapse to the 
real axis as long as the Rayleigh- Taylor is greater than or equal to 0. This allows us to reach 
a regime for which the boundary z develops a vertical tangent. 

Section 5 is devoted to showing the existence of a large class of analytic curves for which 
there exists a point where the tangent vector is vertical and the velocities indicate that the 
curves are going to turn over and reach the unstable regime for a small time. Plugging these 
initial data into a Cauchy-Kowalewski theorem indicates that the analytic curves turn over. 
Therefore the unstable regime is reached. 

Finally, in section 6, a perturbative argument allows us to conclude that we can find 
curves in H'^ close enough to the special class of analytic curves described in Section 5, which 
satisfy the arc-chord and Rayleigh- Taylor conditions. Then we can show the existence of the 
curves passing the critical time and actually turning over. Therefore the unstable regime is 
reached for an entire il^— neighborhood of initial data. 

Remark 1.2 In a forthcoming paper (see J^) we will exhibit a particular initial datum for 
which we will show that once the curve reaches the unstable regime the strip of analyticity 
collapses in finite time and the solution breaks down. In section 8 we provide a very brief 
sketch of our proof of breakdown of smoothness for the Muskat equation. These results were 
announced in f^. 

Remark 1.3 The same approach can be done for the water waves problem, which shows 
that, starting with some initial data given by (a,/o(a)), infinite time the interface reaches a 
regime in which it is no longer a graph. Therefore there exists a time t* where the solution of 
the free boundary problem parametrized by {a,f{a,t)) satisfies \\fa\\L°°{t*) = 00 (see section 
[^. This scenario is known in the literature as wave breaking and there are numerical 
simulations showing this phenomenon 

Remark 1.4 We conjecture that a result analogous to Theorem holds, in which surface 
tension is included. We may simply use the same initial data as in Theorem and take 
the coefficient of surface tension to be very small. The solutions are presumably changed only 
slightly by the surface tension (although we do not have a proof of this plausible assertion). 
Consequently, we believe that Muskat solutions with small surface tension can turn over. 

A similar remark applies to water waves (see theorem \7.1\ ). There exist initial data for 
which water waves with surface tension turn over. A rigorous proof may be easily supplied, 
since local existence (backwards and forward in time ) is known for water waves with surface 
tension (see f^). 

2 The contour equation and numerical simulations 

Here we present the evolution equation in terms of the free boundary which is going to be 
used throughout the paper, and the numerical experiment that motivated the Theorem. 
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2.1 The equation of motion 

By Darcy's law: 

V"*" • u = — (p^ — p^)daZ2{a)S{x — z{a)), 

and Biot-Savart yields 



{P'-P')py f {z{a)-z{a-(3))- 
2tt |z(a) -z(a-/3)|2 



= -'-^^^PV I T~ daz.{a - m. (8) 



using the identity 



For the first coordinate above one finds 

27r U \z{a)-z{a-PW 

27r Jr \z{a) - z{a - I3)\^ 

PV [ a^(ln(|z(a) -z(a-/3)|2))(i/3 = 0. 
Jr 

Therefore 

2vr \z{a) - z{a - P)Y 

Here we point out that in the Biot-Savart law the perpendicular direction appears, but after 
the above integration by parts, we only see the tangential direction. 
Adding the tangential term 

(^2 - p^) f (zi(a)-zi(a-/3)) ^^^ 

PV I 1^^^^ — ^^:^dpdaz{a), 



2vr \z{a) - z{a - /3)| 

we find that the contour equation is given by 

2vr J^\z{a) - z{a- (3)\^ 

For the 27r periodic interface the equation becomes 

z^(a) = ~ r sin(^i(Q) - zi{a - (3)){daz{a) - dazja - 

4vr cosh(z2(a) — Z2{a — (3)) — cos(zi(a) — zi{a — (3)) 

In order to see ([9]) we take z{a) = zi{a) + iz2{a); it is easy to rewrite ([8]) as follows; 

zt{a) = — — PV / — T^dp. 

2-Ki z{a) - z{P) 

The classical identity 

I ^ z \ 1 
7 + E 



-(27r/c)2; 2tan(z/2) 
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allows us to conclude that 

■y (r^\ - ~ f (sinh(z2(Q) - Z2{f3)),-sm{zi{a) - zi{(3))) 
' 47r Jj cosh(z2(a) - z,m - cos(zi(a) - z,m 

where T = R/IttZ. 

Analogously, using the equality 



Att cosh(z2(a) - 2;2(/3)) - cos(zi(a) - 

Att cosh(z2(a) - Z2{p)) - cos{zi{a) - zi{p)) 

and adding the appropriate tangential term, we obtain equation Q. 

2.2 The scenario motivated by the numerics 

Our investigations started with the idea that interesting new phenomena may arise if we 
study three fluids, separated by two interfaces. Careful numerical studies indicated that one 
of the interfaces may turn over. In attempting to prove analytically the turnover indicated 
by the numerics, we discovered that a turnover can occur also for a single interface, i.e., for 
the Muskat problem. This section describes one of our numerical experiments. 

Proceeding as in the preceding section, one can derive the equations modeling the evolu- 
tion of two interfaces separating three fluids with different densities pj (j = 1,2,3). More pre- 
cisely, assume that both interfaces can be parametrized by graphs {a, f{a, t)) and (a, g{a, t)), 
with / lying above g. These equations read in the periodic case, cf. [Ml [15] (this scenario 
has been recently also considered in |20|). 

Ma, t) = pi I[f{;t), f{;t)] + p2l[f{;t),g{;t)], f{a, 0) = /o(a), ^^^^ 

9t{a,t) = p2l[g{-,t),g{-,t)] + pil[g{-,t)j{-,t)], g{a,0) = 5o(a), 

where pj = {pj+i — /3j)/(47r), j = 1,2, and, for given functions u{a), v{a), 

f {dguja) - davja - P)) tan(/3/2)(l - tanh^((^(a) - v{a - /3))/2)) 

l\u,v\ := FV I 7^ dp. 

jf tan2(/3/2)+tanh2((u(a)-z;(a-/3))/2) 

(11) 

The first terms ![/(■, t), /(•, t)] and X[g{-,t),g[-,t)\ in ([TOj) give the velocity of a unique inter- 
face. The cross terms t), t)] and X[g{-,t), /(•, t)] take into account the interaction of 
the two interfaces, and their contribution is getting bigger when the curves are getting closer. 
This, together with the diffusive behavior reported in [16] for the equation 

ft{a,t) = -piZ[f{-,t),f[-,t)], f{a,0) = fo{a), (12) 

and the mean conservation for / and g, motivate the choice of the following initial data, in 
the hope that some non regularizing effect arises from the interaction of the two interfaces; 

q /7r(a - Mi+ri)\ 
0.1 - sin3 (^^ ^ , if a G [Ml - ri,Mi + n], ^^^^ 

0.1, otherwise 



/o(a) 
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2 4 6 2 4 



Figure 1: Left: Solutions to (jlOp with initial data (jl3p -(|14p at times t = (dashed blue), 
t = 3.46 • 10"^ (red points) and t = 7.66 • 10"'' (black). Right: Solutions at t = 1.04 • 10~^ 
(dashed red) and t = 1.84 • 10"^ (black) 

and 

{. 3 /7r(a - Ma + ra) 
sm 
V 2r2 
go{a) = -0.92, 

The choice of parameters Mi = tt + 0.1, ri = 0.7, M2 = 7r/1.2, r2 = 0.3, pi = 207r and 
P2 = tt/20, yielded a strong growth of the derivative in the the lower interface as the two 
curves approach, as shown in Figure [TJ Moreover, after introducing a small modification 
in the lower interface so that the tangent at a certain point becomes actually infinite, and 
evaluating the normal velocity relative to this point along the modified curve, we obtain the 
result plotted in Figure [2j This graphic clearly indicates that the velocity field is forcing the 
interface to turn over. 

The numerical approximation of (jlOp addresses as a main difficulty the absolute lack of 
knowledge about the behavior of the solutions to (jlOp . Indeed, the goal of our experiments is 
precisely the search for some singular behavior. The nonlocal terms make the computations 
expensive and special care has to be taken in order to evaluate the integrands in a neigh- 
borhood of /3 = 0. For this, we used Taylor expansions locally and computed exactly the 
principal value. In this situation, adaptivity is strongly indicated, both in space and time, 
since a good indicator of a singular behavior will be given either by a sudden accumulation 
of spatial nodes or a sudden reduction of the time steps. 

In order to attain the highest resolution in the integration of (llOp and compute the 
solutions shown in Figure [H cubic spline interpolation of the curves f{-,t) and g{-,t) with 
periodic boundary conditions was used. This provides a interpolant of each interface at 
every time and allows, in particular, the evaluation of the convolution terms at any /3 G [0, 27r]. 
Then, adaptive quadrature can be applied to approximate the integrals and evaluate the 
derivative at any time. In the experiments reported, adaptive Lobatto quadrature was used, 
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3 

-0.92, if Q G [M2 -r2,M2 + r2 



(14) 
otherwise. 








2.8 2.9 3 3.1 

Figure 2: Zoom of the interface, modified so that its tangent is vertical at a single point P; 
and the normal velocity along the curve, minus that at P, scaled by a factor of 100. 

by means of the MATLAB routine quadl. For the time integration, the embedded Runge- 
Kutta formula due to Dormand and Prince, DOPRI5(4), was implemented, since the problem 
was not found to be particularly stiff, see for instance [21]. The time stepping was combined 
with a spatial node redistribution after every successful step. For the redistribution of the 
spatial nodes an algorithm following [T7] was implemented, with some modifications taking 
into account that both interfaces are graphs. For several tolerance requirements and different 
choices of the parameters involved in the full adaptive routine, the integration always failed 
at a certain critical time, suggesting the explosion of the derivative at a certain point of the 
lower interface and the lack of validity of (llOp . once this curve stops being a graph. 

The phenomenon described above and the explicit representations of the maximum of the 
solutions derived in [T^, motivated the search for special initial data which allowed us to 
understand that this behavior also arises in the one-interface case. 



3 Instant Analyticity 



Here we show the main estimates that provide local-existence and instant analyticity for a 
single curve that satisfies initially the arc-chord and Rayleigh- Taylor conditions. We consider 
the function 

\z{a) — z{a — /3)p 

and in the periodic setting 



F(z)(a,/3) = — , , , z 7 a,/3GT, 

2(cosh(2:2(a) — Z2[a — p)) — cos(zi(a) — zi[a — p))) 

where ||x|| = dist{x,2Tr'Z). 

If F(z) S then we say that the curve satisfies the arc-chord condition, and the 
norm of F is called the arc-chord constant. 
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Let us clarify the meaning of the above arc-chord condition. Fix t, and assume that 
z{a,t) is a smooth function of a. Suppose F G Letting f3 tend to zero, we conclude 

that \daz{a,t)\ is bounded below. Since also z is smooth, \daz{a,t)\ is also bounded above. 
Consequently, the numerator in the fraction defining F is comparable to the square of the 
arc-length between z{a,t) and z{a — /3,t). On the other hand, the denominator of that 
fraction is comparable to the square of the length of the chord joining z{a,t) to z{a — /3,t). 
Thus, the boundedness of F expresses the standard arc-chord condition for the curve z{-,t) 
together with a lower bound for \daz{a,t)\. 

Theorem 3.1 Let z{a,0) = zo{a) G H'^, F{zo){a,l3) G L°° and daZi{a,0) > (R-T). 
Then there is a solution of the Muskat problem z{a, t) defined for < t < T that continues 
analytically into the strip S{t) = {a + i(^ : |C| < at} for each t. Here, c and T are determined 
by upper bounds of the norm and the arc- chord constant of the initial data and a positive 
lower bound of daZi{a,0) . Moreover, for < t < T , the quantity 

^ / {\z{a ± ict) -{a + ict, 0)^ + \diz{a ± ict)\'^)da 

is bounded by a constant determined by upper bounds for the norm and the arc- chord 
constant of the initial data and a positive lower bound of daZi{a, 0). Above \ -\ is the modulus 
of a complex number or a vector in C^. 

Proof: For the proof we consider the contour z G H'^ with z — (a, 0) periodic and daZi{a, 0) > 
0. In the case of the real line similar arguments hold. The Muskat equation reads 

r sm{zi{a) - zi{a - P)){daz{a) - daz{a - /3)) 
J_^ cosh(z2(a) - Z2{a - p)) - cos(zi(a) - zi{a - p)) 

where we suppose daZi{a,Q) > 0. We also take — = A-n since we are studying the case 
p^ > p^ . For the complex extension one finds 

zJa + iC) = r + ^0 - zija + K- P)){daz{a + Q - daz{a + ij - (3)) 

* J cosh{z2{a -\- iC) - Z2{a -\- iC - P)) - cos{zi{a -\- iC) - zi{a -\- iC - /3)) 

(16) 

We will use energy estimates. Consider 

S{t) = {a + iC G C : Q G T, Id < ct}. 



for c given 



Irll2 

I^I1l2(5) 



{t) = y^ / \z{a ± ict,t) — {a ± ict,0)\^da, 



± 



'^At the end of the proof we can take any c < vavaa(daZ\{a, 0)/|9a2(a, 0)p). 
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where fc > 2 as an integer, and 
F{z){a + iC,f3) 



2(cosh(2:2(a + iC) - ^2(0 + - /?)) - cos(^i(a + «() - zi{a + iC - /3))) ' 

(17) 

with norm 



\\F{z)\\L^iS){t)= sup \F{z){a + iC,P)\. 

a+iCe5(i),/3GT 

Next, we define as follows: 

IklllW = lkllH4(5)W + ll^(^)llL-(5)(i)- 

We shall analyze the evolution of ||-z||_H'4(5)(f). 

Before starting the energy estimates, we mention an idea used previously e.g. in the proof 
of (6.3) in [11]. Suppose A{a,P) is a C^(T) function, and suppose /(a) belongs to L'^{T). To 
estimate 

A{a,a-(3)^cot{^)f{a-f3)dP (18) 
we break up this integral as the sum of 

A{a,a) icot(|)/(«-/3)d/3 (19) 

and 

[[A{a, a -(3)- A{a, a)]^ cot(^)}/(a - (3)d(3. (20) 

The integral in ()19p is simply the Hilbert transform of / and the quantity in curly brackets 
in (llOp is bounded. This idea will be used repeatedly, with A{a, /3) arising from derivatives 
d^z{a,t) up to order 2, and with /(a) = (9^z^(a, t) {fi = 1,2). Whenever we use this scheme, 
we will simply say that "a Hilbert transform arises" . For similar simple ideas used below, we 
refer the reader to the term Ji in pg. 485 in |11) . 

Then, using above scheme, for the low order terms in derivatives, it is easy to find that 

-— / \z{a±ict),t)-{a±ict,0)\'^da<C{\\z\\s{t) + l)''. (21) 
2 dt Jy 

In (j2ip and in several of the estimates, k denotes a enough large universal constant. 
Next, we check that 



1 d 
2dl 



[ \d^z{a±ict,t)\'^da = -^j- / \d^Zj{a ± ict,t)fda 



where 

ld_ 
2di 



I \d^Zj{a±ict,t)\'^da = ^ / d^Zj{a ± ict,t){dt{d'^Zj){a±ict,t) ±icd^Zj{a±ict,t))da. 
T Jt 



22 
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In order to simplify the exposition we write z{a,t) = z{a) for a fixed t, we treat both 
coordinates at the same time, we write (xi, X2) • (2^3, 2^4) = xix^ + X2Xi for Xj € C, j = 1, 4, 
we denote aUct = 7, and we define 

Q{-f,(3) = cosh(z2(7) - ^2(7 - /?)) - cos(zi(7) - zi{j - /3)). 

Then we spht the right hand side of (|22p by writing 

Ii = 5ft / dt4^ ■ dizt{^)da, 

and 



T 



l2 = ^ 94^(7) • icdlz{^)da. 

JT 

In Ji we wiU find the R-T and use it to absorb I2. We wih decompose Ii in order to find the 
terms of at least fourth order. In order to estimate the lower order terms, we refer the reader 
to the paper [TT] (see, e.g.. Lemma 6.1). We have Ii = Ji + J2 + J3+ l.o.t., where 

||1.0.t||^2(T) <C(||z||5 + l)\ 

and Ji, J2, J?, are defined as follows: 



A = ^j^dtz{^) ■ ( y^A(7,/3) ^"^'^^^^(^°^j^^ ^\ d^z{j)-d^z{^-mp)da, 
where ^(7,/?) = cos(zi(7) - 21(7 - /3)), 

J2 = -K ^ dt^) ■ ( I '"'^''|q(^,^'))2~^^^ ^«^(7) - 9.^(7 - /3))i?(7, /3)d/3) d« 
where 

5(7, /3) = (sin(zi(7) - zi(7 - /3)), sinh(z2(7) " ^2(7 " • (^^(7) - ^^(7 - 

and 

We split further Ji = Ki + K2 where 

Ki = '^j^ dUj) ■ (7) [PV ^ ^M|(a,^(^) - a„z(7 - /3))d/3) da, 



^ I dtz{^)-{PV I ^^(d^z{^)-d^z{^-fi))dUi{l-m)da. 



Ail, (3) 

Taking into account the complex extension of the arc-chord condition, it is easy to deal with 
Ki to obtain 

K,<{\\z\\s{t) + 1)'. 
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In K2 it is possible to find a "Hilbert transform" applied to d^zi as in (jlSp . and therefore an 
analogous estimate follows. We are done with Ji. For J2 we obtain similarly 

J2< (11^115(0 + 1)'. 

Next, we split J3 = + where 

We have to be careful, because for real curves is harmless, but for complex curves we need 
to use the dissipative term to cancel out a dangerous term. We denote 

and therefore = Li + L2 where 

Li= / ^{f){^{dizMdlz) + Q{d^z)Q{dlz))da, 
Jt 

L2= [ Q{f){-n{dUMdlz) + Q{dUMdlz))da. 

An easy integration by parts allows us to get 

Li = -l I 'ii{daf)\dtz\^da<C{\\z\\sit) + l)K 

For L2 we find 

L2= ! ^{d^fmdiz)%dlz)da + 2 / %f)%dlzmdlz))da. 

JT JT 

The first term on the right is easy to dominate by C(||z||5 + l)'^. We denote the second one 
by Ml. We claim that 

Ml < C{\\z\\s{t) + 1)'= + KP(/)||H.(5)||A^/'aM|i,(5), (24) 

for K > {) universal constant. To see this, we rewrite 

Ml = -2 / '^{f)%diz)^{K{H{dlz)))da 

JT 

which yields 

Ml = -2 / Ai/2(9(/)9(a^z))K(Ai/2(i/(9^^)))da 
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and therefore 

<C{\\z\\sit) + 1)'= + i^P(/)||H.(5)||AV25M|i.(5). 
Finally we find that 

Ks < C{\\z\\s{t) + 1)'= + K\\'^{f)\\H^s)\\^'^'dU\hisy (25) 

We will use the thickness of the strip to control the unbounded term above. 
For we decompose further: K4 = + L4 + + Lq where 



r^TTT r l^ sin(zi(7)-zi(7-/3)) .xn^5. o^.o. 



dazi{-f) r .2 1 



Inside L3, L4 and L5 we can integrate by parts and therefore 

L3 + L4 + L5 <C(||2||5(t) + l)^ 

In Lq we use the splitting Lg = M2 + M3 where 

M2 = / • 9(A(5^z)) + ^{d^z) ■ ^{A{dtz)))da, 



^3 = - / 5J(-|^)(3f?(a^) • 5?(A(9^z)) + . 9(A(a^)))da. 

In M2 it is easy to find a commutator formula: 

and the appropriate estimate follows. We find that M2 < C(||2;||s + 1)'^. For M3 we write 
M3 = iVi + N2 where 

^1 = - / - mitmidtz) ■ midiz)) + -s{diz) ■ -s{A{diz)))da, 

JT \OaZ\ 

N2 = -m{t)\\A'/\diz)\\l,^sr 
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where 

m{t) = minK(-l^^i^). 
We use the pointwise estimate [lOj 

2gAig) - Kig") > 0. (26) 

Therefore 
as long as 

«(|^) > »■ 

|9«2(7)r 

Remember that mitially jR( jg^^^^^p ) is greater than zero (R-T). We wih prove that it is going 
to keep Hke that for a short time. For I2 we find as before 



I2 = c J^mdtz){j) ■ K(5^z)(7) - K(5^)(7) • Q{dlz){j))da < c\\A'/'{dU)\\h^sy 
Finally 

~ ^ \diz{a ± ict)\'da < C{\\z\\s{t) + 1)^ + {c + K\\Q{f)\\H^s)it) - m(t))||AV2(a4^)|||,(^)(t). 

Note that p(/)||i^2(s)(0) = 0. If c - m(0) < 0, we wih show that 

c+K\\Qif)\\H2^S)it)-mit)<0 

for short time. It yields 

/ \diz{a±ict)\^da<C{\\z\\s{t) + l)'', 
Jj 



1 d 
2dt 



as long as c + i^||Sj(/)||//2(5)(t) — m{t) < 0. Using Sobolev estimates, we proceed as in section 
8 in [11] to show that 

|||F(z)||^o.(s)<C(||z||5(t) + l)^ 



From the two inequalities above and (j2ip it is easy to obtain a priori energy estimates that 
depend upon the negativity of c+K\\Q{f) \\fj2(^s-^{t) — m{t). We get bona fide energy estimates 
as follows. We denote 

MlAt) = Mlit) + l/{m{t) - c - K\\'^{f)\\H^s)it)). 
At this point, it is easy to find that 

-f^MmHHs)it)<C{\\z\\s{t) + l)' 
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using ()23p . and therefore (see section 9 in |TT] for more details) 

J^\\z\\RT{t)<C{\\zUT{t) + lf. 

It follows that 

11.11 ||z||flT(0) + l ^ 

" "'''^^ ^ - (1 - C(|l.bT(0) + 
providing the a priori estimate with C and k universal constants. 
We approximate the problem as follows 

r sin(0, * zl{a) - (t>s * zf{mda{<Pe * z'){a) - d^j^, * ^")(^)) ,o 
^ '^'* J cosh(4(a)-4(/?))-cos(4(a)-4(/3)) ^ 

z^(a,0) = (j)e*ZQ{a), 

where (pei^) = 4' is the heat kernel and e > 0. Picard's theorem yields the existence 

of a solution z^{a,t) in C ([0, r^);iJ^) which is analytic in the whole space for zq satisfying 
the arc-chord condition and e small enough. Using the same techniques we have developed 
above we obtain a bound for z^(a,t) in in the strip S{t) for a small enough T which is 
independent of e. We need arc-chord, R-T, zq £ and c — m(0) < 0. Then we can pass to 
the limit. 



4 Getting all the way to breakdown of Rayleigh- Taylor 

This section is devoted to proving the following theorem. 
Theorem 4.1 Let z(a,0) = z^{a) be an analytic curve in the strip 

S = {a + iC eC : \C\ < h{0)}, 

with h{0) > and satisfying: 

• The arc-chord condition, F(z°)(a + iC,/3) G L°°{S x M) 

• The Rayleigh- Taylor condition, daZi{a) > . 

• The curve z^{a) is real for real a. 

• The functions z^{a) — a and Z2{cy) are periodic with period 2tt. 

• The functions -Zi(a) — a and Z2{oi) belong to H^[dS). 

Then there exist a time T and a solution of the Muskat problem z{a, t) defined for < t <T 
that continues analytically into some complex strip for each fixed t € [0,r]. Here T is either 
a small constant depending only on \\z^\\s or it is the first time a vertical tangent appears, 
whichever occurs first. 

Thus our Muskat solution is analytic as long as daZi{a,t) > 0. 
We will use the following: 
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Lemma 4.2 Let ip{a ± iC) = Ef=-7V ^fce^^°=FfcC. Then, for ( > 0, we have 
d_ 



^ V / \ip{a±iC)\'^da> [ Aip{a ± iC)ip{a ± iQda - 10 [ Aip{a)7p{a)da, (27) 



where A^{a ± iC) = J2k=~N |A:|^fce^'="eT^? . 

Proof: First we shall compute the left hand side in the frequency space: 

N 



[ \^{a±iC)\^da = A7T ^ cosh(2|A;|C). 

4- ■'T 1 »r 



± ■'^ k=-N 

On the other hand we have that 

N 



^ / Aip{a±iC)ip{a±iC)da = 4Tr ^ \k\\Ak\^ cosh{2\k\C) 
^ •^'^ k=-N 



while 

N 



I Atp{a)ip{a) da = 2tt V] |A;||Afcp. 

-'IT U- AT 



k=-N 

Differentiating in ( we obtain 

N 



— / \^{a±iC)\'^ da = 8Tr V sinh(2|A;|C). 

k^N 

The lemma holds since sinh(C) > cosh((^) — 1 for any > 0. 

Corollary 4.3 Let ip{a±i(,t) = J2k=~N ^k{'t)e-^^'^e^^'^ and h{t) > be a decreasing function 
oft. Then 

/ \^{a±ihm^da<^y] I ^^{c^±iHt)Ma±ih{t))da 
-I0h'{t) I Aip{a)Tp{a)da + / ipt{a ± ih{t))ip{a ± ih{t))da. 



^1 _l_ 

This corollary allows us to prove Theorem I4.1[ 

Proof (Theorem 14. ip : The norms ||2||/ffe(5) and \\z\\s are defined as before using the new 
strip S{t) defined by 

S{t) = {a + < eC : Id < h{t)}, 
where h{t) is a positive decreasing function of t. 
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We use the Galerkin approximation of equation (jlSp . i.e. 

9fz[^](C,t) = n^[j[z[^]]](C,t), 

where C £ S{t), IIjv will be specified below, and 

cosh(z2(a) - Z2{l3)) - cos(zi(a) - zi{(3)) 

We impose the initial condition 

z[^](a,0) =z[^l(a). 

Here, for a large enough positive integer N, we define z^^\a,0) from z^{a) by using the 
projection 

oo N 
-oo -N 

We define z^^^{a) by stipulating that 

z[^^(a) — a = n7v[^;i(a) — a] 

and 

zfi(«) = n^[^2°(«)]- 

For large enough, the functions 2^^1(0,0) satisfy the arc-chord and Rayleigh- Taylor con- 
dition. 

We shall consider the evolution of the most singular quantity 



^ / \dtz^^\a±ihN{t),t)\^da, 



where h]\f{t) is a smooth positive decreasing function on t, with /lAr(O) = h{0), which will be 
given below. Also we denote 

SN{t) = {a + iCeC : \C\<hN{t)}. 

^From now on, we will drop the dependency on from z^^l and h]\f{t) in our notation. 
We will return to the previous notation in the discussion below at the end of the section. 
Taking the derivative with respect to t yields 



d f 2 

^ / \dizfj,{a±ih{t),t)\ da 
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23f? / d^Zf,{a±ih{t),t) {dtd^z^{a±ih{t),t) + ih'{t)dlzf,{a±ih{t),t)} da 



/ d^z^{a±ih{t),t) {diUN[J^,[z]]{a ± ih{t),t) + ih' {t)dlz^{a ± ih{t),t)} da 
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/ d^z^{a±ih{t),t) {UN[d^J^[z]]{a ±ih{t),t) + ih'{t)dlz^ia ±ihit),t)} da 



= 2Ky" d^z^{a±ih{t),t) {dtj^[z]{a ± ih{t),t) + ih' {t)d^Zf,{a ± ih{t),t)} da, 

since (9^z^(a =b ih{t),t) is a trigonometric polynomial in the range of Hn- Here = 1, 2. 
Using the above corollary we have that 



d_ 
It 



V f \dizi,{a±ih{t),t)fda<^^S2[A{dU,.){(^±Mt))-dtzt^ia±iHt))da 



T ^ JT 



10/i'(t) / A(a^z/,)(a) ■94z/,(a)da + 2^K / dp^[z]{a,t){a ± ih{t)) ■ d^Zf,{a ± ih{t))da 



We shall study in detail the most singular term in d J[z]{a,t), i.e. 



sm{zi{a±ih{t),t) - zi{P,t)){d^z{a ± ih{t),t) - d^z{^,t)) 



daJ[z]{a±ih{t),t) cosh{z2{a ±ih{t),t) - Z2{/3,t)) - cos{zi{a ±ihit),t) - ziiP,t)/^ 

+ l.o.t = X + l.o.t., 

where ||Lo.t||^2(']p) < + (see [11] and our previous discussion of (fTSl) ). We split 

X in to the following terms 



X = K{a ± ih{t),p){dlz{a ± t) - 9^z(/3, t))d/3 

+ a{a ± t) cot (55^(a ± t) - 5|z(/3, t))d/3 



where 

sin(zi(a,t) - zi{l3,t)) 



K{a,P) 



cosh(2;2(a, i) - Z2{(3,t)) - cos{zi{a,t) - 



daZiia, t) f a — B 

cot ' 



{daZ2{a,t)Y + {daZi{a,t)Y 
and 



a{a, t) 



{daZi{a,t)y + (9a^;2(Q;,t))^' 
Let us denote 

r±(t) = {C€C : C = a± a e T}. 
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Since K{a, j3) is a holomorphic function in a and /3, with o, /3 € S{t), for fixed t we have that 

/■TT 

Xi = / K{a±ih{t),f3)d^z{a±ih{t),t)df3 

J TT 

/■TT 

- / K{a±ih{t),/3)dlz{f3,t)d(3 

J TT 

=^11 + -^^12, 

and integration by parts shows that the term X12 satisfies 1 1-^^121 1 L2(f) <C(||z||s' + l)'^. In 
addition, we can write Xn as follows 



Xn = / K{a ± ih{t),w)d^z{a ± ih{t),t)dw 



'wer±{t) 

I sm{zi{a±ih{t),t) — zi{w,t))d^z{a±ih{t),t) ^ 

Jw&v±{t) cosh(z2(a ± ih{t),t) — Z2{w, t)) — cos(zi(q ± ih{t),t) — zi{w, t)) 

- d^z{a ± ih{t),t)a{a ± ih{t),t)P.V. [ cot f^^±^!M^\ dw 

_py r sm{zi{a±ih{t),t) - zi{P±ih{t),t))d'^z{a±ih{t),t) 

cosh(2;2(a ± ih{t),t) — Z2{(3 ± ih(t),t)) — cos(2;i(a ± ih{t),t) — zi(/3 ± ih{t),t)) 

As before we call 
f{a±ih{t),t) 

^py r s\n{zi{a±ih{t),t) - z^ifi ±ih{t),t)) 

j _^ cosh(z2(a ± ih{t),t) — Z2{P ± ih{t),t)) — cos(zi(a it ih{t),t) — zi{(3 ± ih{t),t)) 

^ p y r sin(zi(a ± ih{t),t) -zi{a± ih{t) - 13, t)) 

cosh(z2(a ± ih(t),t) — Z2{a ± ih{t) — /?, t)) — cos(zi(q; it ih{t),t) — zi{a i: ih{t) — /?, t)) 

Thus 

Xn = d^z{a ± ih{t),t)f{a ± ih{t),t). 
Also we can write X2 in the following way; 

X2 =a{a ± ihit),t) j cot {dlz{a ± ih{t),t) - d^^z{l3, t))d(3 

=a{a ± ih{t),t) [ cot f '^^'Ht)-w \ ^^5^^^ ^ ih{t),t) - dlziw, t))dw 

Jwer±(t) V 2 / 

, . . f / a lb ih(t) — W\ ^r: , 

=a{a ± ih{t), t)P.V. / cot ^-^ dlz{a ± ih{t),t)dw 

JweT±{t) V 2 y 

- a{a ± ih{t),t)P.V. [ cot /^^i^iM^) j dlz{w, t)dw 
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= -a{a±ih{t),t)P.V. j cot (^^^y^) d^zi/S ±ih{t),t)d(3 

= -a{a±ih{t),t)P.V. j ^ csc^ (^^^y^) {d^z{a ± ih{t),t) - d^z{(3 ± ih{t),t))dl3 
and finally 

X2 = - 27ra(a ± t)(A9^z)(a ± ih{t),t). 
Then we find two dangerous terms 



Ii = 2^ I f{a±ih{t),t){d^Zf,){a±ih{t)) ■ {d^z^){a±ih{t))da 

and 



T 



/a = -47r3? / d{a ± ih{t),t)A{d^z^){a ± ih{t)) ■ d^Zf,{a ± ih{t))da. 

The rest can be bounded by C(||2;||5 + l)'^(t) as in the previous section. In order to bound 
Ii and I2 we use the following commutator estimate: 



||A2(/g) - /A2g||i2(Tr) < C||A^+VllL2(T)|bllL2(T), (28) 

for /(a) = Ae^'" and g{a) = E- -N9ke^ '^, where e > and C does not depend on A^. 

The proof of ([28|) will be left to the reader. 

First we estimate /i. We denote 7 = + ih{t). 



Ii=2K / /(7,t)94z^(7,t)a^z^(7,t)da 

5^(/(7)){3?(5^^(7,t))9„(K(5^^(7,t))) + 9(9^^(7, t))9„(9(a^^(7,t)))} da 
-2 / G(/(7)){K(5^^^(7,t))9a(9(5^^(7,t))) + 9(9^,(7, t))9„(5?(5^^(7,t)))} da 



Integrating by parts we have that H-Z^iiHLacir) < C(||z||5 + 1)^. In order to estimate /12 we 
note that /(7,t) is real for real 7. Then 

^{f{a±ih{t),t))=h{t)f±{a,t), 

where 



\\f±\\m{T)<C{\\z\\s{t) + l) 



k 
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Then we can write 

G(/(7))5ft(a^^(7, t))d^{'^{diz^{j, t)))da 

r 

= h{t) r f±{a,m{diz^{j,t))damdiz^{j,t)))da 

J —n 

= -h{t) r f±{a,m{dU^{j,t))AH{Q{dizf,{j,t)))da 

J — TT 

= -hit) r A-2{Ma,midtz^)i^,t))A'2H{Q{dtz^{^,t)))da 

J —n 

= -h{t) r {A^(/±(a,t)»(9^z^)(7,t))-/±(a)A^»(a^^)}A^if(S(9^z^(7,t)))d- 

J — TT 

-h{t) r h{a)Kh^{diz^)K^H{%diz^{-t,t)))da 

J —TT 

<MOI|A^(/±(^i)W^^(.±iMi),i))-/±(^t)A^»(9^Z;.(-±iMi)))llL2(T) 
x\\AhH{^{diz,{- ±ih{t),t)))\\L2(j^ 

+ h{t)\\h\\L^^j-^\\K\^{diz^{-±ih{t)))\\L2^^^^^^ 

Using the estimate (f28]) yields 



a 



<hm^'^'f±\\LHT)\mdiz^{-±Mt)M\LHT)\\^k^^^^^ 

+ Mi)ll/±llL°o(T)l|A^5?('9^M(-±^Mi)))l!L2(T)l|A^9(5^M(-±^Mi))llL2(T) 
< Chm\z\\s + 1)' + Ch{t){Ms + inA'^dXi- ± ih{t),t)\\h(^T) 



= Chm\z\\s + ir + Ch{t){\\z\\s + ir / 94^^(7, AS^^ (7, 

Now Ii is equal to the integral to the left, plus a similar integral that can be bounded in a 
similar way. 

Thus we obtain that 

5^/1 < Ci\\z\\s + 1)'= + Ch{tK\\z\\s + inA'/'dUlhisy (29) 
± 

By assumption the R-T a is bigger than zero for real values. In order to avoid problems 
with the imaginary part we may write 

daZi{a±ih{t),t) daZi{a,t) 

- + h{t)g±[a,t). 



{dazi{a±ih{t))f + {daZ2{a ± ih{t)))^ |9„z(a,t)|2 
where 

\\g±\\miT)<C{\\z\\s + l)''. 
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One finds, 



T \daz{a)\^ 



A{daZ^){a ± ih{t)) ■ d^Zfj,{a ± ih{t))da 



-h{t)2^ / g±{a,t)A{d^Zf,){a±ih{t)) ■ d^z^{a±ih{t))da. 
Jt 

The first term above can be treated as in section 3 taking advantage of the inequality 
26l) . Here we just need daZi{a) > 0. The second term can be treated using the inequality 
as with the term /i. We find that 

5^/2 < Ci\\z\\s + 1)'= + C/.(t)|I<7±|lH2(5)l|Ai/'9^||i,(^), 
± 



and therefore 



^h< C{\\z\\s + 1)" + Chm\4s + l)'\\A'/'dU\l2 



is)- 



(30) 



d 
di 



Using ([29]) and ([30]) we have that 

V / \diz^{a±ihm^da<C{\\z\\s{t) + l)'' -Wh'{t) [ A{diz^){a) ■d^^{a)da 
_i_ Jt Jt 

+{C{\\z\\s{t) + l)^h{t) + ^h'{t)) [ A(5^>)(a ± ih{t)) ■ d^z^ia ± ih{t))da. 
Choosing 

= /i(0)exp(-10C / {\\z\\s + lf{r)dr) 



we eliminate the most dangerous term. The other term in the expression above involves with 
a function on the real line and it is easily controlled. Indeed 

as one sees by examining the Fourier expansion of d^z^{a, t). 
Thus 



10h'{t) / A{diz^){a) ■ d^z^{a)da 
Jt 

And we obtain finally 



< C 



^iiziii<c(iizii5+ir^ 



h{t) 



Jt^l |5^(a ± ihm'da < C{\\z\\s{t) + i; 



fc+2 
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Recovering the dependency on N in our notation we have that 



As in the previous section, we can obtain a bound of the evolution of the arc-chord 
condition that depends on C(||2;[^l + 1)^'+^^ 

This estimate is true whenever t G [0,r/v], where T/v is the maximal time of existence 
of the solution z'^L In addition inequality pip shows that we can extend these solutions in 
H'^{S) up to a small enough time T independent of N and depending on the initial data. 

The above calculation shows that the strip may shrink but does not collapse as long as 
daZi{a, t) > 0. 



5 From an analytic curve in the stable regime to an analytic 
curve in the unstable regime 

In this section we show that there exist some initial data which are analytic curves satisfying 
the arc-chord and R-T conditions such that the solution of the Muskat problem reaches the 
unstable regime. In order to do it we will prove the local existence of solutions for analytic 
initial data without assuming the R-T condition. Then we will construct some suitable initial 
data for our purpose. 

Theorem 5.1 Let zq he an analytic curve satisfying the arc- chord condition. Then there 
exists an analytic solution for the Muskat problem in some interval [— T, T] for a small enough 
T > 0. 

Remark 5.2 Notice that in theorem ^5.1\) there is no assumption on the R-T condition. 
The proof we use here is analogous to the one in 133^ based on Cauchy-Kowalewski theorems 
127\ \28^ (for an application to the Euler equation see f^). Here we cannot parametrize the 
curve as a graph, so we have to change the argument substantially in the proof in order to 
deal with the arc-chord condition. 

Proof: We use the same notation as before. Let {Xr}r>o be a scale of Banach spaces given 
by M^— valued real functions / that can be extended into the complex strip Sr = {a -\- i^ G 
C : \C\ < r} such that the norm 

\\f\\r = y] I \f{a±ir)-ia±ir,0)\^da+ [ \d^f{a±ir)\^da, 

^ JT JT 

is finite and /(a) — (a,0) is 27r— periodic. 

Let z^{a) be a curve satisfying the arc-chord condition and z^{a) G for some tq > 0. 
Then, we will show that there exist a time T > and < r < tq so that there is a unique 
solution to ([IS]) in C([0,r];X^). 



24 



It is easy to check that C Xr' for r' < r due to the fact that ||/||r' < ||/||r- A simple 
appHcation of the Cauchy formula gives 

\\daf\\r'<^\\f\\r, (32) 

for r' < r. Next, we write equation (jl6p as follows: 

zt{a + i(, t) = G{z{a + iQ, t)), 

with 

sin(zi(a + iC)-2i(Q + iC-/3))(9«^(a + 0-5a2;(a + iC-/3)) 
G(z(Q+z4, rjj = / — — — T — —T — — — T^dp. 

cosh(z2(a + zC) - Z2{a + - p)) - cos(zi(a + iC.) - zi[a + iC, — p)) 

We take < r' < r and we introduce the open set O in Sr given by 

= {z,LO(^Xr:\\z\\r<R, \\F{z)\\Loo(^Sr)<R'}^ (33) 

with F{z){a + iC,(3,t) given by (I17p . Then the function G for G . O ^ Xr' is a continuous 
mapping. In addition, there is a constant Gr (depending on R only) such that 

||G(z)||,, <^||z||., (34) 

\\G{z^) - G{z')\\r' < ■^,\\^'' - A\r. (35) 

and 

sup \G{z){a + iO-G{z){a + iQ-P)\<GR\l3\, (36) 

for z,z^ S O. The above inequalities can be proved by estimating as in previous sections. 
Then they yield the proof of theorem 15.11 The argument is analogous to [27] and [28] . We 
have to deal with the arc-chord condition so we will point out the main differences. For initial 
data z^ € satisfying arc-chord, we can find a < < tq and a constant i?o such that 
IIz'^IIj.^ < i?o and 

^ cosh(zO(Q + ip - z^ja + iC- (3)) - cos(z°(a + <) - z?(a + < - /?)) ^ 1 

R-o 



for a + iC, ^ 5,,/ . We take < r < Tq and Rq < R to define the open set O as in (p3 
Therefore we can use the classical method of successive approximations: 



z''+\t) = z^+ I G(z"(s))ds, (38) 
Jo 

for G . O ^ Xr' and < r' < r. We assume by induction that 

||/||,(t) < R, and \\F{z'')\\Lo.^Sr)it) < R 
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for k <n and < t < T with T = mm{TA,TcK) and Tqk the time obtained in the proofs 
in [27] and [28], and Ta determined below. Now, we will check that ||F(z"'"'"^)||^oo(5^)(t) < R 
for suitable Ta- The rest of the proof follows in the same way as in [27], [28j . 
Definitions (j38p and (|17p easily imply that 

|(F(z"+i)(a + iC, 13, t)r'\ > |(F(zO)(a + iC, /?, " ^^.(t' + " if " ^^^^^ + 

To see this, we just use the formulas for cos(a + b) and cosh(a + b), and bounds for the 

r , ■ cosh(a;) — 1 1— cos(a;) sinh(a;) sin(x) r i i i rni r j. i ■ 

functions V — , 3-^, —, — —, tor bounded x. iherefore, takme; 

X X X X JO 

0<r,<n>in{l,^(^-i,)A.}, 

we obtain ||F(z"'''^)||ioo(5^,)(t) < R. This completes the proof of theorem 15. II 

The next step will be the construction of analytic initial data such that 
a. daZi{a) > if a / 0. b. do,zi{Q) = 0. 

c. 5^^2(0) >0. d.daVi{Q)<Q. 

Also zi[a) — a and ^2(0) are 27r— periodic. 

Here Vfj_{a,t), with fi = 1,2, are the velocities given by 

cosh(2;2(a) - Z2[P)) - cos{zi{a) - zi{l3)) 

Notice that in this situation the graph / : M — ?> M defined by the equation Z2{ct) = 
f{zi{a)), has a vertical tangent at the point z{0). See the figure below for an example. We 
shall prove the following lemma: 

Lemma 5.3 There exists a curve z{a) = {zi{a), Z2{a)) with the following properties: 

1. Zi{a)—a and Z2{ot) are analytic 271 — periodic functions and z{a) satisfies the arc-chord 
condition, 

2. z{a) is odd and 

3. dazi{a) >Oifa^O, daZi{0) = and 9^^2(0) > 0, 
such that 



(5„z;i)(0) 



5„ 



sin(zi(a) - zi{/3)) 



^ cosh(2;2(a) - Z2{l3)) - cos{zi{a) - zi{l3)) 



{daZi{a) - daZi{(3))df3 



< 0. 



a=0 

(39) 
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-4-3-2-1 1 2 3 4 



Proof: We shall assume that z{a) is a smooth curve satisfying the properties [2] and [3l 
Differentiating the expression for the horizontal component of the velocity, it is easy to obtain 

cosh(z2(a) - Z2{a - p)) - cos(zi(q) - zi(a - p)) 

cos(zi(a) - 2i(a - P)){daZi{a) - dgZija - j3)f 
cosh(2:2(a) — Z2{a — j3)) — cos(zi(a) — zi{a — /?)) 

^ r sin(zi(a) - zi{a - - dlz^{a - 

cosh(z2(a) — Z2{a — (3)) — cos(zi(a) — zi{a — /?)) 



sin((zi(Q) - zi{a - (3))){daZi{a) - daZi{a - /?)) 
sinh(z2(a) - 2:2(0 - l3)){daZ2{a) - daZ2{a - 



(cosh(22(a) — -22(0: — /?)) — cos(zi(a) — zi{a — 

sin((zi(Q) - zi{a - f3))){daZi{a) - daZi{a - /?)) 
sin(z2(a) - Z2{a - l3)){daZi{a) - daZi{a - /?)) 



dp 



(cosh(z2(a) — 2^2(0: — — cos(zi(a) — zi{a — /S)))^ 
Evaluating this expression at a = we have that 



dp. 



(r. ^, , cos(zi(/3))(a^zi(/3))^ + sin(zi(/3))agzi(/3) ^^ 

cosh(z2(p)) - cos(zi(/3)) 
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■ ( (RWP) ^^^ sm(zi(/3))(9aZi(/3) - sinh(z2(/?)) (^0^^2(0) -daZ2{(3)) 

Sm[Zi[p ) )OaZl{p ) . . , , f n\\\'2 

(cosh(z2(/3)) - cos{zi{f3))y 



Integration by parts yields 

sin(zi(/3))a2zi(/3) 



^ cosh(z2(/3)) - cos(zi(/3)) 



dl3 



cos zi /? — . . r~n^'^P 

cosh.{Z2{p)) - cos(zi(/?)) 



+ / sm(zi(/3))5c,Zi(/3) taw /fl^^^2 '^f^' 



, sin(zi (/3))a^zi(/3) + sinh(z2(/3))a aZ2 (/?) 
(cosh(z2(/3))-cos(zi(/3)))2 

The above integrals converge because zi and Z2 satisfy the properties 2 and 3. Therefore 
we obtain that 



(0) /' sin(..(^.inh(..(m ^^^^ ^ ^^ 

(cosh(2;2(/3)) - cos(zi(/3)))^ 



o;) f sin(zi(^))smh(2;2(/3)) ^ ,^ ,.„s 

= 29^22(0 / 7 daZi{P)dp (40 

Jo (cosh(z2(/3)) - cos(zi(/3)))^ 

^From the expression (j40p we can control the sign of {daVi){0). In order to clarify the 
proof we shall take 

zi(/3) = -sin(/3) + /3. 

We construct the function 22 in the following way: 

Let /?! and /32 be real numbers satisfying < /3i < /32 < vr, and let z*{P) be a smooth 
function on [— 7r,7r], with the following properties, 

a. is odd. b. {d[sz*){0) > 0. 

c. z*{f3) > if /3 € (0,/3i). d. z*{/3) < if /3 G (/3i,/32] 

e. <Oif^e [/32,7r]. 

For a positive real number b to be fixed later, we define a piecewise smooth function z{(3) on 
[— vr, vr], by setting 

m = bz*i/3) if|/3|</3i, 
m = z*ip) if /3i < |/3| < vr. 

Then 



(cosh(z(/3)) -cos(zi(/3)))2 
is negative and independent of b, while 

^■^1 sin(zi(/3))sinh(5(/3)) 
(cosh(5(/3))-cos(zi(/3)))^ 



/ 

JO 



5a^i(/3)d/3 



28 



tends to zero as 6 — oo. 

Therefore, we can fix b large enough so that 



sin(.,(/3))smh(i(/?)) .^^^^^ ^ ^ 



(cosh(z(/3))-cos(zi(/3)))2 

It is now easy to approximate z{/3) in L^[— vr, vr] by an odd, real-analytic 27r— periodic function 

r sin(.i(/3))sinh(z2(/?)) 
Jo (cosh(z2(/3)) - cos(zi(/3)))^ 

and (9aZ2(0) > 0. 

The conclusions of lemma (j5.3p follow, thanks to ([I0|) . 

Theorem (15.111 and lemma (15.31) allow us to show the breakdown of the R-T condition. 



Theorem 5.4 Let zq a curve satisfying the requirements of lemma l[5.3\) . Then there exists 
an analytic solution of the Muskat problem satisfying the arc-chord condition in some interval 
[— T, r] such that for small enough T > we have that: 

1. daZi{a, —t) > Va and 

2. daZi{0,t) < 

for all t G (0, r]. In addition daZ2{0,t) > in [—T,T]. 

Proof: We use theorem (j5.ip to obtain the existence and from lemma (15. Sp we have that 

{dtdaZi)iO,0) <0. 

Remark 5.5 For t G [— T, 0], our solution satisfies 

nimdaZi{a,t) > c\t\. 

a 

This follows easily, since daZi{a,0) has a non- degenerate minimum ata = 0, and dtdo,zi{0,0) < 
0. 

6 From a curve in in the stable regime to an analytic curve 
in the unstable regime 

Finally we show that there exists an open set of initial data in the topology satisfying 
the arc-chord and R-T conditions such that the solution for the Muskat problem reaches the 
unstable regime. This section is devoted to proving theorem (II. ip . 

Proof of theorem (jl.ip : The idea is simply to take a small i/^-neighborhood of the initial 
data of an analytic solution. Let zq be a curve as in lemma (15. Sp . Let z{a, t) with t € [—T, T] 
for some T > 0, the solution for the equation 1^ given by theorem 15. 11 We consider the curve 



29 



w^{a) = {{wg{a)i, (i(;|(a))2) which is a small perturbation in H^{T) of the curve z{a,t) at 
time t = —6, with < 5 < T, i. e. 

\\wU-)-z{.,-6)\\H. = M\\m<e. 

Also, Wg{a) satisfies the R-T condition 

(^wiia) = {p^ - p'^)da{wl)i{a) > 0, 

if < (5 < Jo and < e < £{S). From now on, we take e and 5 to satisfy this condition. Also, 
we may take £{6) < Eq. 

Since z{a,0) = zo{a) is a smooth curve satisfying the arc-chord condition we can assume 
that there exist £q > and < 5o < T such that 

sup \\wl{-)\\Hi(T) <C{zo,£o,6o), (41) 

0<£<£o,0<5<5o 

and 

sup ||F('u;f)||i^(Tr) < C(zo,eo,5o). (42) 

0<£<£o,0<<5<<5o 

Now, let the curve w^{a,t) be the solution to the equation 

sm{{w')i{a,t)-{w')i{^,t)) 



cosh(('u;=)2(Q,t) - {w^)2i(3,t)) - cosiiw')iia,t) - (w;^)i(/3, t)) 
X {da'w^{a,t) - da'w'^{/3,t))d/3 
w^{a, —5) = tt^f (a). 

^From theorems 13.11 and 14.11 and the inequalities (j4ip and ()42p we see that we can choose Eq 
and 6 small enough in such a way that w^{a, t) is well defined, for all < e < Eq, in t S [—6, 0] 
unless w^{a,t) loses the R-T condition. That means there exist some point ao and some time 
to S [—5,0] satisfying aw^iao,to) < 0. Also, for small enough eo and fixed 5 we have that 

CTw^ (a) >a> 0, (43) 

where a is a real number independent of £. The numbers eo and 6 are fixed for the rest of 
the proof. 

If there exist times t such that there exists some point ao with a"^£(ao,t) = 0, we denote 
the first of these times to be G (— 5, oo). Also we set T'^ = min{r^,0} and P = [—5,T^]. 
Due to (f43l) we have that 

inf >h> -6, 

0<e<e{S) 

for some number tb. 

^From the proof of theorems 13.11 and 14. 1 1 we know that there exists a function h{t), given 
by the expression 

Ca{t + 5) -5<t<ta 
Ca{ta + d)e~^'^^'-'^^ t>ta, 



hit) 
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where ta (small enough), Cq and Ca are constants which only depend on the constant C{zQ,eQ, 6) 
(see ()4ip and (I42p ). such that w^{a,t) is an analytic function in the strip 

5(t) = {CGC : m)\<h{t)}, 

and also 

{\\w^.,t)\\s + lf<Ca, 

for some large enough k and t G [ta>^^] (notice that the constants ta, Ca and Ca do not 
depend on e). 

In this situation we claim the following: 

^ / \d^{w'{a±ih{t),t) - z{a±ih{t),t))\'^da 

<C{\\dt{w'{-±ih{t),t)-zi-±ihit),t))\\l,^^^ + \\w'{. + ih{t),t)-zi- + i^^^^ (44) 

for t ^ P and where C is a constant just depending on C{zq, eq, 5). 

We will prove this inequality at the end of the section. Let us assume that ()44p holds. 

We notice that we can always choose either a subsequence {en}i° with e„ — > when 
n — )■ CO such that T*"" < Vn or a subsequence {em}i° with — >■ when m — )■ oo such that 
2^em > Q \/jy^ {ih.e case in which there exist only a finite number of times can be treated 
as this last case). We deal with these two cases, I and II, separately: 

I. < for all e. From inequality (|44p we can take e small enough such that 



has norm < Ce in H'^{S{T')). 
Note that 



= mmdc{w')i{a,T'') > -Ce + mindaZ^ > -Ce + c\T'\ 

a a 

by the remark at the end of section 5. Thus, \T^\ < Ce. 
Then 

z{a, T^) — ZQ{a) 
has norm < Ce in H'^{S{Q)) and therefore 

\{daMi){ao,T') - {d^{v,,)im\ < Ce, 

and we can conclude that 

{dc.Mi){ao,T') <0. 

Here we recall that 

sin((w;^)i(a,t) - {w^)i{P,t)){do,w^{a,t) - daW^{P,t))dp 



cosh(('u;=)2(a, t) - {w^)2{li, t)) - cos{{w^)i{a, t) - {w^)i{l3, t)) 



Applying the same argument as in section [5] to the curve w^{a,T^) we finish the proof of 
theorem 11.11 in the case T"^ < for all e. 
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II. > for all e. Then we can apply a Caucliy-Kowalewski theorem to the initial data 

■w'^{a, 0) — z{a, 0) 

satisfying 



- z\\s(o) < Ce. 



For t > small enough, z{a, t) is in the unstable regime. We achieve the conclusion of 
theorem 11.11 by continuity with respect to the initial data. 

The rest of the section is devoted to proving inequality We shall denote 7 = a + ih{t) 
and d{-~f,t) = d'^{w{'^,t) — z{'y,t)) (we omit the superscript e in the notation) and we recall 
that w{a, t) and z{a, t) are real for real a (therefore we obtain similar similar estimates for 
7 = Q — ih(t)). In order to prove inequality (jH]) we have to compute the following quantity 

d 



j \d{-^,t)\^da ='^^y d{-f,t)dt{j,t)daj +2^^ih'{t) J d{j,t)dad{-f,t)da 

Again we treat in detail the most singular term in dt{'y,t). Recall K{a,(3) from section 3 and 
write Kyj and for corresponding expressions arising from z and w. Then we have that 

dt(7, t) = r KUl, 7 - ^)dUwil, t) - w{j - /3, t)dp 

K,(7, 7 - f3)dl{z{^, t) - z{j - (3, t))dp + l.o.t(a, t), 

J —IT 

where 



/ d{-f,t)l.o.t{a)da'> < C(||d(- + i/i(t), f)||22(Tr) + I k(- + ^/i(t), - ^(- + «ft(t), t)| li2(T))- 



Here C is a constant which just depends on Eq and 6. 
We can write 



K^{-f, 7 - (3)dl{wi-f, t) - w{-f - /?, t))d(3 

K,(7, 7 - P)dU4l, t) - z{j - (3, t))d(3 



+ / {K,(7,7-/?)-if^(7,7-/3)}9^(z(7,t)-z(7-/3,t))d/3 



K^j, 7 - /3)5„(d(7, t) - (i(7 - /3, t))d/3 
+ r {K,{j, 7 - /3) - K^il, 7 - /3)} dlizil, t) - z(7 - /?, t))d/3 

J —TT 

Xi{a,t) + X2ia,t). 
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Therefore 

d_ 
di 



23?<! / d(^Xi{a,t)da\ +2^{ I d{j,t)X2{a,t)da 



+ 2^<ih'{t) / d{-f,t)dad{-f,t)da 



Following the computations in section 3 when t G [—5, ta] and those in section 4 when t S 
[ta,T^] we have that 



I \d{j,t)fda < C\\d{- + ih{t),t)\\l,^^^ + 



(i(7, t)X2(a, t)da 



In addition 



sin(tf;-^(7) — ^^(7 — 



cosh(ti;2(^) — tt;2(-y — /3)) — cos(w^(7) — ^1)^(7 — /?)) 
sin(zi(7)-^n7-/3)) 



cosh(z2(^) _ ^2^^^ _ ^-j-) _ cos(z^(7) — ^^-(7 — /?)) 

(5„«;i(7))2 + (5««;i(7))2''° V2 



i^«.(7,7-/3) 
i^.(7,7-/3)- 



9a^^(7) 



cot 



/3 



+ 



cot 



/3 



(5„u;i(7))' + (5a^n7))2 (5.^1(7))' + (5a^H7))=^ J V2 



cot 



/3 



Also, 



tan ( — 
2 



\dlz{a ± ih{t),t) - dlz{a ± ih{t) - /3,t)\ < C 
since z is the analytic unperturbed solution. Therefore 

We are done. 

7 Turning water waves 

Let us consider an incompressible irrotational flow satisfying the Euler equations 

p{vt + v-Vv) = -Vp-gp{{),l), (45) 
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where p satisfies (j2l3p and = 0. This system of equations provides the motion of the 
interface for the water wave problem (see [3l [25] and references therein), whose contour 
equation is given by 

zt{a^ t) = BR{z, ijj){a, t) + c{a, t)daz{a, t), (46) 

and 

iOt{a,t) = -2dtBR{z,uj){a,t) • daz{a,t) - da{^^Q-^){a,t) + da{cuj){a,t) ^^^^ 
+ 2c{a, t)daBR{z, uj){a, t) ■ daz{a, t) — 2gdaZ2{a, t). 

The values of z{a, t) and 10(0, t) are given at an initial time tQ-. z{a, Iq) = z^{a) and w{a, to) = 
w^{a). For more details see |12| . 

As an application of section 5, we can consider initial data given by a graph (a,/o(a)) 
and show that in finite time the interface evolution reaches a regime where the contour only 
can be parametrized as z{a,t) = {zi{a,t), Z2{a,t)), for a € M, with daZi{a,t) < for a € /, 
a non-empty interval. This implies that there exists a time t* where the solution of the free 
boundary problem reparametrized by {a,f{a,t)) satisfies ||/a||L°° (i*) = 00. 

Theorem 7.1 There exists a non-empty open set of initial data z^(a) = (a,/o(a)) and 
w^{a), with /o € and G H'^, such that in finite time t* the solution of the water wave 
problem ^J^^^J^ given by {a, f{a,t)) satisfies ||/a||L°° (^*) = 00. The solution can be continued 
for t > t* as z{a,t) with daZi{a,t) < for a £ I, a non-empty interval. 

Proof: Let us consider a curve z*[a) G satisfying 1., 2. and 3. of Lemma 15.31 We 
point out that analyticity is not required here. In order to find a velocity with property 
(I39p we pick for water waves u}{a,t*) = —daZ2{a) and a suitable z{a,t*) = z*{a) as an 
initial datum. Notice that the tangential term does not affect the evolution. Then, with the 
appropriate c{a,t), we can apply the local existence result in [12]: There exists a solution of 
the water wave problem with z{a,t) € C{{t* - 6,t* + 6];H^), uj{a,t) € C{[t* - 5,t* + 6];H^) 
and (5 > small enough. The initial data promised by theorem 17.11 are any sufficiently small 
perturbations of z(a,t) and w{oL,t) at time t = t* — 6. 

8 Breakdown of Smoothness 

In [5] we will exhibit a solution z{a, t) of the Muskat equation, with the following properties. 

1. At time to, the interface is real-analytic and satisfies the arc-chord and Rayleigh- Taylor 
conditions. 

2. At time ti > tQ, the interface turns over. 

3. At time ^2 > ^i) the interface no longer belongs to C^, although it is real-analytic for 
all times t € [toj ^2)- 

In this section we provide a brief sketch of our proof of the existence of such a Muskat solution. 

Our Muskat solution z{a,t) will be a small perturbation of a Muskat solution z^^(a,t), 
with the following properties. 
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4. z^^{a,t) is real analytic in a, for |9a| < e^^ and |r| < r'^''. 

5. For t G [— r^'^jO), z^^{a,t) satisfies the Rayleigh- Taylor and arc-chord conditions. 

6. For t = 0, the curve z^^(a,t) has a vertical tangent at a = 0. 

7. For t € (0,T^'^], the curve z^^{a,t) fails to satisfy the Rayleigh- Taylor condition. 

This paper constructs Muskat solutions z^^ satisfying 4., 5., 6. and 7. Our problem is to pass 
from z'^^ to a nearby Muskat solution z satisfying 1., 2. and 3. The idea is as follows. 

So far, we have studied the analytic continuation of Muskat solutions to a time- varying 
strip 

S{t) = {\Qa\ < h{t)}, 

in the complex plane. In our forthcoming paper [5], we will study the analytic continuation 
of a Muskat solution to a carefully chosen time- varying domain of the form 

n{t) = {\'^a\ <h{?fia,t)}, (48) 

defined for t G Here, r is a small enough positive number. 

For t G [—t^^,t], we will work with the space H^{^l{t)), consisting of all analytic functions 
F : Q{t) 1-^ C'^ whose derivatives up to order 4 belong to L'^[dQ,{t)). 

We will pick our time-varying domain in (I48p so that h{x,t) > for all {x,t) G 
M/27rZ X [-Ti°,r) and h{x,T) > for ah x G M/27rZ \ {0}, but /i(0,r) = 0. Thus, the 
domain Q{t) has 'thickness' zero at the origin. Consequently, H'^{Q,{t)) is not contained in 
C^(R/27rZ). 

We will also take r < t^^ and h{x,t) < e^°, so that the Muskat solution z^^{a, t) continues 
analytically to i^{t), for each t G 

We can therefore pick an 'initial' curve z^{a), such that 

8. z^{a) — z^^{a,T) belongs to H'^{^1{t)) and has small norm, yet 

9. z°{a) does not belong to C"^(R/27rZ). 

We solve the Muskat problem backwards in time, with the 'initial' condition 

10. z{a,T) = z°(a). 

By a more elaborate version of the analytic continuation arguments used in this paper, we 
find that our Muskat solution exists and continues analytically into for all t G [t*,T] 

(for a suitable time t^,); moreover, 

11. z{a,t) — z^^{a,t) has small norm in H^{Q{t)), for all t G [t*,T]. 
Here, either 

12. = -r^o or 

13. a modified Rayleigh- Taylor condition, adapted to the time- varying domain, fails at time 

We can rule out 13., thanks to 11., together with our understanding of z^^{t) and Q{t). 

Thus, we obtain a Muskat solution z{a,t), satisfying 9., 10., 11. and 12. Properties 1., 2. 
and 3. of z(a, t) now follow easily. 
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